ULB-TH/07-04 



Bound on the Dark Matter Density in the Solar System from 

Planetary Motions 

J.-M. Frere, F.-S. Ling, G. Vertongen 

Service de Physique Theorique 
Universite Libre de Bruxelles 
B-1050 Brussels, Belgium 

PACS numbers : 95.10.Ce, 95.35.+d 

Abstract 

High precision planet orbital data extracted from direct observation, spacecraft explorations and 
laser ranging techniques enable to put a strong constraint on the maximal dark matter density of 
a spherical halo centered around the Sun. The maximal density at Earth's location is of the order 
10 5 GeV/cm 3 and shows only a mild dependence on the slope of the halo profile, taken between and 
-2. This bound is somewhat better than that obtained from the perihelion precession limits. 

Introduction 

The existence of Dark Matter in the Universe has been inferred from the observation of structures 
much larger than the Solar System. Since Fritz Zwicky first suggested in 1933 the presence of a non- 
luminous matter component in the Coma galaxy cluster pQ, numerous experimental observations as 
well as theoretical studies have fortified this paradigm. On the cosmological scale, the analysis of the 
cosmic microwave background mapped by the WMAP satellite [2] has led to the so-called Concordance 
Model of cosmology j3j, in which only around 4% of the total energy density is made of ordinary matter, 
while the rest is shared between Dark Matter (26%) and Dark Energy (70%). The need for Dark Matter 
has also been shown in structure formation scenarios and simulations. On the galaxy cluster scale, the 
amount of Dark Matter inferred from galaxy kinematical observations matches the strong lensing data 
of the Hubble Space Telescope [4]. In the next future, strong and weak lensing surveys with SNAP and 
LSST should provide us with a detailed sky map of the gravitational shear induced by Dark Matter [5] . 
On a galactic scale, numerous spiral galaxies have been studied and yield flat or even rising rotation 
curves, that extend well beyond the visible disk [5] . This is usually interpreted as due to a Dark Matter 
halo, much bigger than the disk itself. 

If Dark Matter is an essential component of the universe, there is some hope to discover it in our 
own galaxy, the Milky Way. As expected, modeling and understanding the detailed structure of the 
halo becomes crucial to estimate the intensity and the spectrum of the Dark Matter signal in direct 
as well as indirect detection experiments. It is generally accepted that the local Dark Matter density 
at Earth's location should be around 0.3 GeV/cm 3 , although a much higher value is possible in the 
presence of clumps, or near a caustic [7j- Therefore, one can imagine that the Solar System itself is 
surrounded by a local subhalo, with a local density well above the galactic value. Several processes 
have been invoked to clump the Dark Matter in the Solar System. For example, Dark Matter particles 
could become trapped inside the Sun's gravitational potential through multiple scatterings [8] or they 
could be captured in the Solar System as a result of the gravitational pull from planets [§] . It might 
even be that the existence of the Solar System itself is evidence for a local subhalo. In any case, a 
Dark Matter halo centered around the Sun will influence the motion of the planets. The purpose of 
this paper is to put limits on such a possibility. 

It is worth emphasizing that the presence of Dark Matter is NOT needed on the scale of the Solar 
System, despite all the aforementioned evidences at larger scales. Celestial Mechanics is one of the 
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greatest achievement of newtonian mechanics, that enabled to calculate ephemerides tables with high 
precision, predict eclipses, and even led to the discovery of Neptune in 1846 and of Pluto in 1930. 
Of course, we now know that newtonian mechanics is not the whole story, as people looked for a 
hypothetical planet to explain Mercury's perihelion precession, but never found it... 

Nowadays, planet orbits have been so precisely measured [10| that they are able to constrain the 
local Dark Matter density. Previous estimates give bounds of the order of 10 -19 to 10~ 20 g/cm 3 , and 
are derived from observational limits on the secular perihelion precession of the inner planets (i.e. the 
residual precession of a planet after all known effects, like the general relativity contribution, or the 
precession due to the influence of other planets are subtracted) ( [12] , [13] , [14] ) . In this paper, we point 
out that the presence of a Dark Matter halo centered around the Sun will not only cause the planets 
to precess, but will also alter all the keplerian parameters used to fit a planet's orbit. In particular, 
if the halo density is high enough, the third Kepler's law cannot hold for all planets, even when all 
the uncertainties of each planet's orbital parameters are taken into account. Current planetary orbit 
determinations are precise enough to provide a meaningful constraint. Of course, at such precision, a 
planet orbit cannot be viewed as a pure keplerian elliptic orbit anymore, as the mutual influence of 
planets introduces many periodic or secular corrections to the keplerian parameters. For our purpose, 
it is however enlightening to continue to discuss the problem in terms of keplerian parameters, in the 
spirit of a perturbation expansion, where the different perturbations are added separately. 

Hereafter, we will therefore describe the orbit of a planet, which is approximately elliptic, by only 
three parameters, which are the semi-major axis a, the eccentricity e, and the revolution period T . In 
the presence of a Dark Matter halo, we first notice that the third Kepler's law can always hold for 
one planet, irrespective of the density of the halo (but still small enough compared to the mass of the 
Sun so that perturbation scheme is valid), because the orbital motion can only probe the total mass 
interior to the orbit. However, with several planets, the halo density is directly limited by the precision 
at which the third Kepler's law is verified. 

The paper is organized as follows. In section 1, we summarize the salient features of the keplerian 
orbit. We recall the various relations between the dynamical variables and the measurable parameters. 
In section 2, we show how the dynamical variables get modified by the introduction of a radial per- 
turbation. In particular, we calculate the shift in the gravitational constant AK and the shift of the 
perihelion AO. Our method also provides a useful formula for the perihelion precession. In section 3, 
these analytical results are applied to a heliocentric spherical Dark Matter halo and numerical estimates 
are derived and discussed. 



1 The keplerian orbit revisited 

The relative motion of a planet of mass mp around the Sun (with mass Mq) is usefully parameterized 
by 6 parameters: K = G(Ai© + mp) gives the strength of the keplerian potential Uo(r) = —K /r, 
Lq is the angular momentum per unit mass, Eq is the energy per unit mass (in absolute value), the 
semi- major axis a and the eccentricity e characterize the elliptical orbit, and T is the revolution period 
around this orbit. All of them are not independent, as 

A'., = (1) 

* = §, P) 
L% = K p, (3) 

where p = a(l — e 2 ) is called the semilatus rectum. The subscript under a parameter is used to 
designate its value in the unperturbed scheme (i.e. without halo). For given values of Kq, Lq, and Eq, 
the orbital elements a, e and the period T can be calculated. The period T is found by computing 

T = 2 dr , (4) 

Jr min y/-2(Eo + U Q (r))-L*/r* 
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where r min = a(X — e) and r max = a(l + e) are respectively the orbit's perihelion and aphelion distances 
to the Sun. In this integral, r min and r max are calculated as the roots of the function 

f (r) = -2(E + U (r))-^ . (5) 

The third Kepler's law (eq. (HJ) is obtained by using the following parameterization of the elliptical 
orbit 

r = a(l — ecos£) , (6) 

where £ is called the eccentric anomaly parameter. 
To calculate the angle variation around the orbit, 

e = 2 f m " Lo/r2 dr (7) 

° Jr min ^/-2(E + U (r))-L*/r* ' 

the parameterization 

r= l + fcostf ' (8) 

where 8 is the polar angle, also called the true anomaly parameter, proves to be more convenient. For 
a pure keplerian motion, Qq is of course 2tt, as there is no perihelion precession. 
The planets parameters are summarized in table [TJ 

2 Shift of the keplerian parameters due to a perturbation in 
the potential 

In this section, we will derive how the keplerian parameters are affected by a perturbation in the 
potential, 

U(r) = U (r) + U 1 (r) , (9) 

which is supposed to be small, U\(r) <C Uo(r) along the unperturbed orbit. 

As discussed in the introduction, the orbital elements a, e and the period T of the planets in the 
Solar System are known to a high precision. We will estimate in section 3 the uncertainty on our results 
allowed by a small variation of these parameters. In this section, however, we will keep a, e and T 
fixed. As a consequence, the dynamical parameters E, L and K will necessarily undergo a shift after 
the introduction of the perturbation. We will write 

E = E Q + AE , 
L 2 = L 2 + AL 2 , 
K = K Q + AK, 

and 

f(r) = -2{E+ U{r)) - ^ = / (r) + A/(r) . (10) 

Notice that AE and AL 2 can be expressed in terms of U\(r) and AK by requiring the invariance of 
the orbital elements. The period is then calculated as a function of AK, 

T(AK) 9 dr T A f^ AK ^ dr nil 

T(AK) = 2 L. vmrmm ~ L, • 

where the last integral is converging by construction. By requiring the invariance of the period, one 
finds 

AK= — [{l + e) 2 U 1 {r max )-{l-e) 2 U 1 {r mm )}+ / V?(r(£)) cos £ d£ , (12) 

2e irae J Q 

with V%(r) = r 5 U\(r), and the derivation is taken with respect to r. 
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Our methodology also enables us to calculate the perihelion shift per revolution, 



ae K ^ _ r~ . 2 r ulm} r i ie)mS m . (13) 

It is worth noticing that the perihelion shift is independent of AK, i.e. independent of a variation 
in the keplerian potential, as expected. Furthermore, eq. (|13p is equivalent to the Landau & Lifschitz 
formula [15] 

Ae = 7TT (t- f ' r 2 (6)UM0)) de) , (14) 



8Lq \Lq j 

although this may not be obvious at first sight, which we checked numerically. 



3 Constraints on a heliocentric Dark Matter halo 

In this section, we examine the consequences on the planetary motions of the presence of a spherically 
symmetric Dark Matter halo centered around the Sun, with a density profile defined by 

p(r) = Po (^) 7 , (15) 

where we take < 7 < 2, and ro is the mean Earth to Sun distance, so that po is the Dark Matter 
density at Earth's location. 7 = corresponds to a constant density halo. Note that such halo will 
affect the orbits of the planets only if it is heliocentric. 7=1 corresponds to the NFW profile that 
appears systematically in Dark Matter simulations [15] . 

The halo mass enclosed in a sphere of radius r is equal to 



M(r)= 4?rP0 _ 7 . (Hi) 

(3 - 7Fo 

It will correspond to a gravitational potential 



AnGpo 
(3-7)(2-7K" 7 



CMO = - ; ^ ■ (17) 



The shifts of AK and the perihelion given by eqs. ([T2^) and (fT3")l reduce to 

Aif(po,7) = -(4- 7 )GM(a), (18) 
Ae( Po ,7) = -tt(3-7)^. (19) 

Notice that eq.(19) differs from the results for 7 = in [H],[T3], but agrees with [14] . 

These shifts provide a constraint on the Dark Matter density at Earth's location when compared to 
the experimental limits found in table [21 Notice that the values listed for AO and Aa are the result of 
a general x 2 orbit fitting of a huge amount of astronomical data, and therefore cannot be considered as 
independent from each other. It is curious that the residual perihelion shift for Venus is not compatible 
with zero. Also, there is no data on Neptune's perihelion precession due to the fact that Neptune's 
revolution period is quite long compared to the time span its orbit has been observed, i.e. 90 years. 
The rationale for deriving a constraint from the calculated values of AK goes as follows. Essentially, 
the maximal shift is limited by the precision at which the period and the semi-major axis of a planet 
are known, by application of the third Kepler's law. Now, the relative precision of these two elements 
are reasonably of the same order of magnitude, as they arise from the same fitting procedure. Also 
notice that the masses of the planets are determined with great accuracy by spacecraft probes [17] . 
Therefore, the maximal shift is limited by the orbit precision combined with a possible modification 
of the solar mass. As already mentioned in the introduction, the knowledge of the orbital elements a 
and T for one planet is not sufficient to constraint the halo density, as the shift AK can always be 
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Table 1: Planets keplerian parameters a, e, T and Kp = Gmp (from [17]). Notice that GM & = 
1.327 10 20 m 3 s- 2 . 



Planet 


Period (days) 


a (AU) 


e 


K P (m^s-' 2 ) 


Mercury 


87.96935 


0.38709893 


0.20563069 


2.20 10 i3 


Venus 


224.70096 


0.72333199 


0.00677323 


3.25 10 i4 


Earth 


365.25696 


1.0000001124 


0.01671022 


4.04 10 i4 


Mars 


686.9601 


1.52366231 


0.09341233 


4.28 10 13 


Jupiter 


4335.3545 


5.20336301 


0.04839266 


1.27 10 17 


Saturn 


10757.7365 


9.53707032 


0.05415060 


3.79 10 w 


Uranus 


30708.16002 


19.19126393 


0.04716771 


5.79 10 15 


Neptune 


60224.9036 


30.06896348 


0.00858587 


6.84 10 15 



Table 2: Experimental limits (from [TO] . [TT] , [T7] ) 



Planet 


AG (arcsec/century) 


Aa (m) 


AK P (m d s- 2 ) 


Mercury 


-0.0036(50) 


0.105 


9.14 10 s 


Venus 


0.53(30) 


0.329 


6.0 10° 


Earth 


-0.0002(4) 


0.146 


2.45 10 Y 


Mars 


0.0001(5) 


0.657 


2.8 10 5 


Jupiter 


0.0062(360) 


639 


9.68 10 s 


Saturn 


-0.92(2.9) 


4222 


1.08 10 9 


Uranus 


0.57(13.) 


38484 


7.59 10 9 


Neptune 


no data 


3463309 


1.4 10 1U 



absorbed in a modification of Mq. For several planets, however, the difference of the measured value 
of if is a direct probe of the halo mass. 

In figure 1, the maximal allowed Dark Matter density at Earth's location from the K shift constraint 
and from the perihelion bounds is plotted against 7. We can see that the first constraint is somewhat 
more stringent. We also notice that the maximal density does not change much when 7 is varied 
between and 2, because it is determined by the precision of the Earth's and Mars' orbits, which are 
the most constraining. 

Finally, let us discuss the so-called Pioneer anomaly in the context of a heliocentric Dark Matter 
halo. It consists of an unexplained constant and uniform acceleration 

a P = (8.74 ± 1.33) ■ 1(T 10 m/s 2 . (20) 

directed towards the Sun and detected by spacecrafts Pioneer 10 and 11 at a distance between 20 
AU and respectively 70 and 40 AU [T5] from the Sun. Such an acceleration would correspond to a 
halo profile characterized by 7 = 1 and po ss 8 • 10 9 GeV/cm 3 , which would have dramatically altered 
the planetary orbits in the solar system. The Pioneer anomaly cannot therefore be explained with a 
spherical halo with monotonically decreasing density. 

Conclusion 

We have evaluated the bound on the maximal allowed Dark Matter density at Earth's location from 
planetary data using the perihelion and the K shifts. The latter bound is somewhat better, although 
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Figure 1: Maximal allowed Dark Matter density at Earth's location, as a function of the halo profile 
parameter 7 



such precision has to be considered as quite optimistic. Conservative bounds could be one order of 
magnitude higher. We have also derived a useful expression for the computation of the perihelion shift. 
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